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Abstract. We study the problem of counting the total number of affine so- 
lutions of a system of n binomials in n variables over an algebraically closed 
field of characteristic zero. We show that we may decide in polynomial time 
if that number is finite. We give a combinatorial formula for computing the 
total number of affine solutions (with or without multiplicity) from which we 
deduce that this counting problem is #P-complete. We discuss special cases 
in which this formula may be computed in polynomial time; in particular, this 
is true for generic exponent vectors. 



1. Introduction 

A binomial ideal in the ring k[xi, . . . , Xn] of polynomials with coefficients in 
a field k, is an ideal generated by binomials: ax"' — bx^, where a,/? £ N" and 
a,b € k* . Binomial ideals are quite ubiquitous in very different contexts particu- 
larly those involving toric geometry and its applications ^1 1^ , in the study of 
semigroup algebras, and in the modern versions of hypergeometric systems of dif- 
ferential equations |25l [7| . While binomial ideals are quite amenable to Grobner 
and standard bases techniques ^1 EHI , they also provide some of the "worst-case" 
examples in computational algebra, such as the Mayr-Meyer ideals |22| . 

In this paper we consider ideals generated by n binomials in R := k[xi, . . . , x„], 
with char(fc) = 0. Let k denote the algebraic closure of k. We are interested in 
determining when the number of solutions in fc" is finite and non zero (i.e., when 
the given binomials define a complete intersection in R) and, in this case, to count 
the number of solutions, with or without multiplicity. We will obtain properties 
of these ideals directly in terms of the given data: the exponents a,/3, and the 
coefficients a, b. 

Our starting point is then a system of n binomials in R, with non-zero coefficients. 
Thus, we may assume that they are of the form 

(1.1) pj(c; x) := x""^ - Cjxf^^ ; j = 1, . . . , n, 

where aj,(3j G N", aj ^ (3j. Let J be the ideal generated by pi, . . . ,p„ in the 
polynomial ring fc(c)[x]. Given a choice of coefficients c G (fc*)", let Jc be the ideal 
in R generated by pi(c; x), . . . ,p„(c, x) and Vc C fc" the variety defined by jTc- 

Proposition l2.1l which is a restatement of results in ^U], gives a complete picture 
of the number of solutions of the system Hl.l|l in the algebraic torus (A;*)". Let B 
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\ OLn- fin j 

whose j-th row is the vector aj — (3j. Then, for generic coefficients c G (fc*)", 
Vc n (fc*)" consists of | det-Bj-many points all of which have multiplicity one (this 
may be seen directly or as a simple instance of Bernstein's theorem). In fact, if 
deti? ^ 0, this is true for all c £ (fc*)". On the other hand, if dct B = 0, then, 
for coefficients c £ (fc*)" not satisfying the algebraic conditions H2.2|l it holds that 
Vc n (fc*)" = , while if the coefficients satisfy (j221l, the variety ("1 (fc*)" has 
codimension equal to the rank of B. We set 6 := | det B\. 

Deciding whether the system (|1.I|) has a non-empty, finite set of solutions in fc" 
is more involved. We must, first of all, consider the possibility that some exponent 
vector aj or /3j may vanish. This is equivalent to the statement that some variables 
Xj are invertible modulo the ideal J. The reduction to the case when this does 
not happen is accomplished in Proposition 12.51 We may then assume that £ Vc 
for all choice of coefficients. Now, in the generic case det B ^ 0, Theorem 12.61 
gives a condition on the exponents of the system that guarantees that the system 
(11.1(1 is a complete intersection for all c £ (fc*)". If, on the other hand, dct_B = 0, 
Theorem 12 . fil only implies that ((1.1(1 is a complete intersection for a generic set of 
coefficients c £ (fc*)". Indeed, in this case, algebraic conditions such as ((2.2(1 enter 
into play. This leads to the notion of generic complete intersection, that we will 
abbreviate by gci. We will say that pi , . . . , p„ is a gci if J7c is a complete intersection 
in R, i.e., Vc is a finite non empty set, for generic coefficients c £ (fc*)". 

Even though Theorem l2 . 61 gives a combinatorial criterion for deciding if pi , . . . , p„ 
is a gci, its verification requires 2" steps. One of the main results of this paper is 
Theorem 12.121 where we describe a polynomial-time algorithm to decide whether 
Pi, . . . ,Pn is a gci directly from the exponents aj , (3j . 

Given a generic complete intersection pi, . . . ,p„, let 

(1.3) d := dinifc fc[a;i, . . . ,x„]/Jc ; D := dimk k[xi, . . . ,Xn]/ \/jc 

be the total number of points in the variety Vc, counted with and without multi- 
plicity. Given an index set L C {1, . . . , n}, we denote by /xl, the number of points 
in V(J) n fc£, kl := {a; £ fc" : = if and only if £ e L} , counted with 
multiplicity. We set [n] := {1, . . . , n} and fi := the multiplicity at the origin. 

In Section 121 we compute d, D, and fiL for a gci. A key ingredient is what we 
call parametric reduction, which allows us to reduce the study of generic complete 
intersection binomial ideals to a particular class of ideals with a normalized presen- 
tation. We show in Theorem 13 . 21 that we can keep track of the various multiplicities 
through the process of parametric reduction. We then compute d and D for so- 
called irreducible systems. We show that an irreducible system that is in normal 
form may behave in one of three possible ways: its binomials are a standard ba- 
sis for either a global or a local term order, or they are weighted homogeneous. 
This allows us to read off the dimension and multiplicities from the exponents (cf. 
Theorem 13.5(1 . Interestingly, the linear algebra problem that underlies these re- 
sults appeared in the work of Vinberg about Cartan matrices JHl Theorem 4.3]. 
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For generic exponents, a binomial system in normal form is irreducible and has 
det i? ^ 0. Hence, Theorem 13.51 gives a polynomial time algorithm for computing 
the number of solutions of a complete intersection binomial system with generic 
exponents and arbitrary non-zero coefficients. 

We next consider the case of a general gci. Using a well-known quadratic-time 
algorithm, due to Tarjan '30', we find a block decomposition of the system into 
irreducible ones. From this decomposition we construct an acyclic directed graph 
naturally attached to the system. In Theorem l3 . 1 51 we give an explicit combinatorial 
formula to compute the dimensions and multiplicities of the system from this graph. 

Section ^ is devoted to counting complexity issues. We reverse the correspon- 
dence from binomial systems to acyclic digraphs and assign to each such graph a 
simple binomial system. The number of solutions of this system corresponds to 
invariants of the graph whose computation is known to be #P-complete. Indeed, 
we show that particular instances correspond to counting independent sets in bi- 
partite graphs, or more generally, antichains in a poset; both of these problems are 
known to be #P-complete ,3 It t24j . Hence, even though the problem of deciding 
whether a system is a gci as well as the problem of counting the number zeros in 
the torus of the binomial system defined by are solvable in polynomial time, 

we prove in Theorem 14 . 31 that counting the total number of affine solutions, with or 
without multiplicity, is a #P-complete problem. Thus, binomial systems furnish a 
very simple example of the type of problems, "easy" to decide but "hard" to count 
that motivated Valiant's introduction of the notion of counting complexity 13 Ij . Fi- 
nally, in Proposition 14.51 we identify another class of systems whose solutions may 
be computed in polynomial time. 

The last section of the paper is devoted to a brief discussions of some of the 
applications of this work which motivated our study. We show, first of all, how 
Theorem 13. 151 may be applied to compute the multiplicity and geometric degree 
of the primary components of a lattice basis ideal J C fc[a;i, . . . , Xm]- This, in turn, 
may be used to describe the holonomic rank of Horn systems of hypergeometric 
partial differential equations and to study sparse discriminants, generalizing the 
codimension-two case. 8 ,7. Finally we recall the results of 29, Chapter 10] relat- 
ing the study of systems of partial differential equations with constant coefficients 
with that of the corresponding algebraic system. 

2. Complete Intersections and normal eorms 

We begin by considering the question of when binomials pi{c; x) , . . . , Pn (c, x) 
as in ()1.1() define a complete intersection when viewed as elements of the Laurent 
polynomial ring S :— k[xf^, . . . ,xf^^]. Let B be the n x n exponent matrix defined 
in H1.2|l . We note that even though the rows of B are only defined up to sign, this 
will not affect our arguments. It follows from ^JJJ Theorem 2.1] that if deti? ^ 
then, for any choice of coefficients in (fc*)", pi{c;x), . . . ,pn{c',x) define a regular 
sequence in S. Moreover, the system of equations 

(2.1) Pj{c;x) = ; j = l,...,n 

has I det-Bj-many solutions in the algebraic torus (fc*)" and all of them are simple. 

On the other hand, if deti? = then pi(c; x), . . . ,p„(c; cc) does not define a 
complete intersection in S for any choice of coefficients. Indeed, if the system (I2.1|l 
has a solution x G (fc*)", it will necessarily have infinitely many. Let TZ be the 
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lattice of relations 

n 

n := {to e Z" : ^ mj{aj - (3j) = 0}. 

For any to G 7?. we have a fc*-action on the set of solutions of (|2.1|) defined by 
{t; x) 1-^ {f"''^xi, . . . , t™"^Xn), and therefore the set of solutions could never be finite. 
Note also that if deti? = then, for generic coefficients Cj, (|2.1|l has no solutions. 
In fact, ii X € (fc*)" is a solution of (|2.1|) we have 

^.aj-ft ^ g^.^ foraU j = l,...,n, 

and therefore 

n 

Yl c"/' = 1 , for all men. 

Thus, a 1/^ , . . . ,1/'^ is a basis of TZ, a necessary condition for pi(c; x), . . . ,p„(c, x) to 
have a solution in (fc*)" is that, 

" i 

(2.2) Yi ^7 = 1 foi' all £=l,...,r. 

3 = i 

This condition is also sufficient. Suppose that H2.2|l holds and let C be the sub- 
lattice of Z" spanned by aj — (3j, j = 1, . . . ,n. Denote by p: £ — > fc* the group 
homomorphism (i.e., the partial character) defined by 

The equalities in (|2.2|) imply that p is well-defined and, since up to a monomial 
(which is invertible in the Laurent polynomial ring), 

it follows from pITl Theorem 2.6] that pi{c;x), . . . ,p„(c, x) define an ideal in S of 
codimension equal to the rank of C. Hence we obtain: 

Proposition 2.1. Let pi(c]x), . . . ,pn{c,x) he as in and B as above. For 

any choice of coefficients c G (fc*)", the ideal they generated in S is a complete 
intersection if and only if det B ^ 0. If det B = and the identities 12. ^) are 
satisfied then the binomials 1^1.1]) define an ideal in S of codimension equal to the 
rank of B. 

In the remaining part of this section, we will discuss criteria for deciding when 
Pi, . . . ,pn is a gci. Since we are not assuming that supp(aj) nsupp(/?j) = 0, where, 
for V e R": 

supp(w) := {i e [n] : ^ 0}, 

the matrix B, by itself, does not allow us to recover the exponents of the binomials 
(11.1(1 . It is useful to introduce the following concept, already present in the work 
of Scheja, Scheja, and Storch I2(jj : 

Definition 2.2. Let pj = x"^ — Cj x^^ , j — 1, . . . ,n, he a system of binomials in 
k[xi, . . . ,Xn]. For each index set KG [n], let 

(2.3) Z{K) {j e [n] : supp(aj) n K 7^ and supp(/3j) n K ^ 0} . 
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We start by showing that we can restrict ourselves to the case where G Vc. 
Since this property is equivalent to the statement that all exponent vectors are non 
zero, it is independent of the choice of coefhcients. We want to identify all indices i 
for which Xi is invertible modulo the ideal ^J, i.e., the Xi coordinate of any solution 
to the system of binomials is necessarily non zero. Set /q = and, for ^ > 1, let 

li ■= IJ'f™PP("j) ■ '5^PP(^i) ^ ^^-i} U'L^^PP(^j) ■ supp(q;j) C li-i} 



and / = li . Induction on £ shows easily that if i S /, the variable Xi is invert- 
ible modulo the ideal and, conversely, that these are all the variables invertible 
modulo J. Thus, after reordering of variables and polynomials, we may assume 
that the variables Xr+i, . . . ,Xn are invertible and that the binomials Ps+i, . . . ,Pn 
involve only the variables Xr+i, . . . ,a;„, while for j < s both monomials x"^ and 
x^^ are divisible by at least one of the variables Xi, i < r, i.e., that ^([t"]) = [s]. 
Following ^2 we define: 

Definition 2.3. Let x' :— {xi, . . . , Xr), c' := (ci, . . . , Cr). For j < s, set 
(2.4) Pj{c';x') =pj{c';{xi,...,xr,l,---,l))- 

Then, the binomial system {pi, . . . ,ps} C k{c')[x'] is called the derived system of 
pi, . . . ,Pn. We denote by 13 the associated s x r matrix as in p.2|l . 

Note that G V(pi, . . . ,ps) and that the matrix B is of the form 

= ( 4 

Lemma 2.4. Assume pi, . . . ,p„ as in is a gci and let r, s be as above. Then, 

r = s and det(i32) ^ 0. 

Proof. Since the variables Xr+i, . . . ,Xn are all invertible modulo J7, the system of 
equations ps+i = • • • = = 0, is equivalent to the system x'^^^^^ — Cj, for all j = 
s + 1, . . . , 71. Hence, arguing as in the discussion leading to Proposition 12. II we see 
that each integer relation among the vectors aj — (3j,j = s + 1, . . . ,n imposes a 
polynomial condition on the coefficients as in (|2.2|) . If s < r, then n—r < n—s and so 
there exists a non trivial relation. Therefore, pi, . . . ,p„ has generically no solutions, 
a contradiction. On the other hand, if s > r, or if r = s and det(_B2) = 0, then, 
generically, the system ps+i(2^r+i, ■ ■ • , Xn) — ■ ■ ■ = p„(a;,.+i, . . . , n) =0 has either 
no solutions or infinitely many in (fc* )""'". Since any solution of these equations 
may be extended to a solution of (12.111 by setting xi = ■ ■ ■ — Xr = 0, we get a 
contradiction again. So s = r and dot (52) ^ 0, as claimed. □ 

Proposition 2.5. Letpi, . . . ,pn, B be as above. Assume that s = r and det(_B2) 7^ 
0. Let pi, . . . ,pr be the derived system. Then pi, . . . ,p„ is a gci if and only if 
pi, . . . ,pr is a gci. 

Proof. Assume pi,...,p„ is a gci and let U be an open dense subset of (fc*)" 
such that the binomials with coefficients in U define a complete intersection ideal 
in k[xi, . . . ,Xn]- It suffices to show that the intersection of U with the fiber 
(k*Y X {(1, ...,!)} is also Zariski dense in the fiber. Let a" G (fc*)""'' be such 
that U n {{k*Y X {a"}) is Zariski dense. Let A" G (A:*)""'' be a common zero of 
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Pr+i{a"; x), . . . ,pn{a"] x). Then, since s = r, the change of variables that sends Xi 
to itself for i = 1, . . . ,r and 

Xj^Xj/X'j, j = r + l,...,n, 

transforms any of the last n — r polynomials pj, j = r + 1, . . . ,n, into a non-zero 
multiple of — x^^ and, for i < r, the binomial Pi into a non-zero multiple of 

where a'l,(3'l € N""'' denote the vectors consisting of the last n — r coordinates 
of UijPi. Since this scalar transformation in the coefhcient space {k*Y preserves 
Zariski dense subsets our assertion follows. 

Conversely, assume that pi,...,Pr is a gci and that det(i32) ^ 0. Let (p be 
a non zero polynomial such that y(c') ^ for a given r-tuple of coefficients 
c' = (ci, . . . , Cr) implies that the corresponding polynomials pi{c';x'), . . . ,Pr{c'; x') 
define a complete intersection. Denote as before c" = (c^+i, . . . ,c„) and consider 
the rational function 

If tjj{c',c") is defined and non zero, then for any choice of the | det(i32)|-inany roots 
A" of the last n — r polynomials, the specialized system 

Pi(c';(x',A")) = ---=p.(c';(x',A"))=0 

has finitely many solutions and, consequently, pi, . . . ,p„ is a gci. □ 

The following result is a reformulation of Theorem 2.3 in 

Theorem 2.6. Let pi, . . . ,pn be as in and suppose that £ Then, 

pi, . . . ,pn ia a gci if and only if \Z(K)\ < \K\ for all K C [n]. 

Proof. Suppose there exists K C [n] such that |Z(iir)| > \K\. Assume that K is 
maximal with this property. After reordering, if necessary, we may assume that 
K = {r +l,...,n} and Z{K) = {s + 1, . . . ,n} where s < r. Since G V( J), 
the maximality assumption implies that the first s binomials depend only on x' ~ 
(xi, . . . ,Xr). Otherwise, we may assume that there exists ki > r, fci G supp(Q;s). 
Since G V(J'), supp(/3s) ^ 0. If there exists ^2 > r, fc2 G supp(/3s), then s G Z{K) 
which is a contradiction. Therefore, supp(/3s) C [r] and for any £ G supp(/3s), 
K' := K U {1} satisfies Z{K) U {s} C Z{K'). Hence \Z[K')\ > \K'\ and this 
contradicts the maximality of K. 

Thus, for a given choice of coefficients, the system 

(2.5) p^{c;x')^---^ps{c;x')^0 

is either inconsistent or its solution space has dimension at least r — s > 0. Since, 
any solution of H2.5|l can be extended to a solution of the full system by setting the 
if-coordinates equal to zero, it follows that pi, . . . ,p„ is not a gci. 

Conversely, suppose |Z(iir)| < \K\ for all K C [n]. In order to show that 
pi, . . . ,Pn is a gci it suffices to prove that given any subset L C [n] , for generic 
coefficients pi{c;x), . . . ,pn{c,x) has at most finitely many solutions with zeros in 
where 

(2.6) kl = {xek"^ -.xe^Oii and only if ^ G L} . 
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Assume that for some choice of coefficients, there exists a sohition in fc£. Then, for 
any i ^ Z{L), pi depends only on the variables in J, the complement of L in [n] 
and hence, since G V(i7), Z{LY C Z{J). Since, by assumption |^(i)| < \L\ and 
\Z{J)\ < \J\, we deduce that 



and therefore = \L\. Reordering we may assume that J = Z(J) = [r] and 

let Bi{L) denote the r x r exponent matrix as in (|1.2|l . If deti3i(L) = 0, then 
for generic coefficients the first r binomials have no solutions in {k*Y and hence, 
generically, pi{c;x), . . . ,pn{c,x) have no solutions in On the other hand, if 
deti?i(L) ^ then, for all choices of coefficients in (A;*)'', there exists finitely 
many solutions of pi = • • • = p,- = in {k*Y and hence finitely many solutions of 
Pi(c;a;), . . . ,p„(c, x) with zeros exactly in L. □ 

Remark 2.7. Note that in the proof of Theorem 12.51 we have shown that if 
Pi,...,p„ is a gci, L C [n], and fc£ is as in 1)2. 6|l . then, for generic coefficients, 
there exists a solution in if and only if |2'(L)| = \L\ and, after reordering so 
that Z{L) = L {r + 1, . . . , n}, the binomials pi, . . . ,pr depend only on the first 
r variables, and the corresponding r x r exponent matrix Bi{L) is non-singular. 
Moreover, for generic c G (fc*)", there are | det i?i(L)j-many points (counted with- 
out multiplicity) in Vc n fc£. Then, the number of points in Vc, counted without 
multiplicity, is given by 



where /i/, is the total number of points in Vc H fc£ counted with multiplicity. We will 
develop in Section|3|the combinatorics needed to describe all sets L with /i^ / and 
we shall show in Section 0] that counting the number of such sets is a #P-complete 
problem. 

Note that if G V( J') , the condition that pi, . . . ,pn is a gci depends only on 
the combinatorics of the exponents aj , Pj . It follows from Proposition 12.11 and 
Theorem 12.61 than, when det(-B) 7^ 0, if pi, . . . ,p„ is a gci, then it is a complete 
intersection for any choice of the coefficients (as long as Cj G fc*). 

The variant of the Fischer-Shapiro criterion embodied in Theorem l2 . fil allows us to 
determine whether pi, . . . ,p„ is a gci. However, this involves checking exponentially 
many conditions, one for each subset K C [n\. We will now show how this can be 
done in a number of steps that depends polynomially (on n). We begin with the 
following simple corollary to Theorem 12. 61 

Corollary 2.8. Suppose pi, . . . ,p„ is a gci and G "V{J'). Let 

M = ; j = l,...,7i} 

denote the set of monomials appearing in pi, . . . ,pn. Then for each i G [n] there 
exists ri > such that x^' G M. . 

Proof. If for some i G [n], x\^ ^ M for all > 0, then Z{{1, . . . , z, . . . , n}) — [n], 
contradicting Theorem 12.61 □ 

One can easily give examples showing that the necessary condition in Corol- 
lary is not sufficient to guarantee that pi, . . . ,p„ define a gci. However, the 
following stronger notion provides a sufficient condition. 



\L\ < \Z{jr\ < \Z{L)\ < \L\, 



(2.7) 




Mi 7^0 
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Definition 2.9. We say that pi, . . . ,p„ are in normal form if and only if for all 

i G [n] 

Pi = - Cix'^' ; ri > 0, (3i ^ 0. 
Note that if the system is in normal form then G V(J'). 

Proposition 2.10. Assume pi, . . . ,p„ are in normal form. Then pi, . . . ,p„ is a 
gci. 

Proof. For any K C [n], Z{K) C K and the result follows from Theorem 12. 61 □ 

We will next show how to reduce ourselves to systems pi, . . . ,p„ in normal form. 

2.1. Parametric Reduction. Let pi, . . . ,pn be a binomial system and suppose 
that they satisfy the necessary condition in Corollarv l2.8l but that it is not possible 
to relabel variables and binomials, or invert the coefficient of one or more binomials, 
so as to put the system in normal form. This means that one of the binomials must 
contain two monomials of the form and x^j^ with i ^ j. Then, after relabeling 
we may assume that p„ is of the form 

(2.8) p„ = - c„a;;^_i , l,m>0. 

Let q :— gcd(m,£) and set m' := m/q, £' := £/q. We will consider the polyno- 
mial map that sends polynomials in n variables xi, . . . ,Xn to polynomials in n — 1 
variables mi, . . . , Un-i- 

(2.9) Xi^Ut, i^l,...,n-2; Xn-i ^ u^^-i] a:„ i-^ 

Let pj, i = 1, . . . , n — 1 be the image of the binomials pi, . . . ,p„_i. We will refer 
to pi, . . . ,Pn-i as a parametric reduction of pi, . . . and denote by J the ideal 
they generate in k{ci, . . . ,c„„i)[ui, . . . ,u„-i]. 

Proposition 2.11. Suppose pi, . . . ,Pn-i is a parametric reduction of pi, . . . ,pn 

and let B and B be the associated matrices. Then \ det-B| — q ■ \ det-Bj. Moreover, 
G V(J^) if and only if d V( »/) and, in this case, pi, . . . ,pn is a gci if and only 
ifpi, . . . ,p„-i is gci. 

Proof. The matrix B is of the form 

_ bl ... bn-2 bn^l bn 

\ ... £ J 

where 5i, . . . , 5„ are vectors in Z"^^. On the other hand, the matrix B is given by 

B= ^bi ... bn-2 i'K-l + m'bn j 

The first assertion now follows from a last-row expansion of det B. 

Suppose now that pi, . . . ,p„ is not a gci. By Theorem 12.61 there exists K C [n] 
such that \Z{K)\ > \K\. U K C [n - 1], then Z{K) C [n ~ I] as well and therefore 
by Theorem 12. 61 Pi . . . . ,p„_i is not a gci either. If n G if, then taking K — K\{n} 
we get that Z{K)\{n} C Z{K). Hence |^(^)| > \K\ and pi, . . . ,Pn-i is not a gci. 

Conversely, if pi, . . . ,p„_i is not a gci then there exists K C [n — 1] such that 
\Z{K)\ > \k\. If X C [n - 2] we take K = K and then Z{K) = Z{K); if, on the 
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other hand, n — 1 G K, then we take K = KU{n} in which case Z{K) = Z{K)\j{n}. 
In either case \Z{K)\ > \K\ and we are done. □ 

The resuhs of this section may be summarized in a polynomial-time algorithm 
to check whether a binomial system is a gci. 

Theorem 2.12. We may decide in polynomial time whether pi^ ... is a gci. 
Moreover, if it is known that det _B ^ we can check if pi{c] x), . . . ,p„(c, x) is a 
complete intersection in time 0(n^). 

Proof. It is easy to see from the procedure for constructing the derived system 
that this step may be accomplished in at most 0{n^) steps. If the number of non- 
invertible variables does not equal the number of binomials in the derived system 
then, by Lemma |2. 41 pi, . . . ,p„ is not a gci. Again by Lemma |2. 41 we next check 
whether det ^2 7^ (this is, of course, unnecessary if it is known that det B ^ 0). If 
so. Proposition 12 . 51 allows us to restrict ourselves to the derived system. We move 
down the list of binomials searching for binomials of the form — cx^^' . Whenever 
such a binomial is found wc do parametric reduction and reduce by one the number 
of binomials and of variables. This step is then repeated until there are no longer 
any binomials of that form. Clearly, this process stops after a quadratic number of 
steps. Then pi, . . . ,Pn is a gci if and only if CoroUarv 12 . 81 holds . This verification 
can certainly be carried out in quadratically many steps. □ 

Example 2.13. Consider the following binomials in k[xi, . . . ,xs]: 

Pi^xj-xl; p2 ^ X1X2 - X1X3; P3 = xlx2X3 - xl; 

P4 = xl~xlxl; p5 = a;^-xg; pa ^ XsXq - X2X3XjXs; 

P7 = X5XT - Xj; ps ^ xl ~ xixexrxs, 

where, since deti? ^ 0, we have set all coefficients Cj — 1. Although the system 
satisfies the necessary condition in Corollary 12.81 it is not in normal form. We 
may apply parametric reduction simultaneously to the binomials pi and p^ by 
considering the polynomial map from k[xi, . . . , xg] to k[ui, . . . , ug] that sends: 

xi^ul; X2^u\] X3 ^ U2; X4^U3; 

X5 ^ u\] Xq i-^ U4; Xj 1-^ M5; Xs ^ Uq. 

Here we have taken into account that the gcd of the exponents in p^ is 2. After 
changing signs when necessary, the new system pi, . . . ,pq is in normal form: 

Pi^ul-ulu2; P2^ul-u\u2\ 

P5 = M5 - W4U5; Pq=uI- uluiU^UQ. 

Thus, we conclude that pi, . . . ,ps defines a complete intersection. We will compute 
the numerical invariants of this system in Example 13.171 

3. Computing the number of solutions 

We recall that if pi, . . . is a gci then we denote by d (respectively D) the 
number of points in Vc H fc" counted with multiplicity (respectively without multi- 
plicity), for a generic choice of non-zero coefficients. Similarly, recall that for any 
index set L C {1, . . . , n} we denote by /i^ the number of points in Vc H k^ counted 
with multiplicity, where fc£ is the set of points in affine space whose coordinate 
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xe = precisely when £ d L. In particular, /i — /i[„] denotes the multiplicity at the 
origin. 

If Pi, . . . ,p„ is a gci but ^ V(J'), then it follows from Lemma [2.41 and Propo- 
sition [^21 that the invariants d and D of pi, . . . ,p„ are obtained from those of the 
derived system by multiplying times | det B2 \ ■ We will assume from now on that 
no variable is invertible modulo J', i.e., that G V(i/). 

We begin this section by showing that it is enough to compute the desired nu- 
merical invariants d,D,^L, for ideals in normal form. We then show that if the 
system is irreducible, in a sense made precise below, then the only zero outside the 
torus is the origin and its multiplicity may be easily computed from the exponents 
of the system. Finally, we consider the general case and show how the various 
dimensions depend on the combinatorics of the irreducible components. 

3.1. Multiplicities and parametric reduction. Suppose pi(c; x) , . . . , Pn{c, x) is 
as in with p„ — xfj — c„x™_]^, £,m > 0. Let q = gcd{£, m) and 

n' — T^' — r' T™' 

We will denote by d' , D' , ^'^ the corresponding invariants for pi, . . . ,p„_i,p^. 

We show, first of all, that by keeping track of q we may assume without loss of 
generality that m and t are coprime. 

Lemma 3.1. With notation as above, set m' = m/ q, I' ~ ijq, = a;^ — c^a;™_]^ 
and let B and B' be the corresponding matrices. 

(1) |detS| = q- |detB'|. 

(2) pi, . . . ,p„ is a gci if and only i/pi, . . . ,p„_i,pj, is a gci. 

(3) For any index set L d {I, . . . ,n}, fiL = 'Z'Ml- 

(4) d^q-d' and D = q- D' . 

Proof. The first assertion is trival while the second one follows from Theorem 12.61 
In order to prove assertion 3, let (ci, . . . , c„) € (/c*)" be such that Jc is a complete 
intersection and decompose 

(3.1) Pn — xl^ — CnX^_i — (x^ — S,x'^_i) , 

where Wq denotes the q-th roots of c„. For any A S Vc, we have 
dinife {Rx/{Jc)x) - ^i^-t (^a/(J5)a) , 

where := (pi(c;x), . . . ,p„_i(c; x), xf^ - ^x'^Li). Therefore, 

dimfc {Rx/{Jc)x) = Y {Rx/iJOx) , 

By a scalar change of variables it follows that 

Y dim,. {Rx/{Ji)x) , 
A6V(j5)nfc2 

is independent of ^ G Wq and, since it agrees with /i^, we obtain that 

IJ-L = qfJ-L 
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as claimed. The last assertion follows directly from the previous one and the fac- 
torization □ 

We next show that multiplicities are not altered under parametric reduction. If 
the binomial system pi, ... ,pn is a gci, and p„ — xf^ — Cnx'^_i, £,m > coprime, let 
Pi, . . . ,Pn-i be the binomial system obtained through parametric reduction. We 
will denote by d, D and /i^ the corresponding invariants. 

Given L C [n] we denote by L := L fl [n — 1]. Conversely, given L C [n — 1] set 
L = Lifn — l^L and L = LU {n} otherwise. Note that if L C [n] is such that 
/iL ^ then either L C [n — 2] or both n — l,n G L. Hence, the correspondence 
L ^ L establishes a bijection between index sets L C [n] such that fiL ^ and 
subsets L C [n — 1] such that /i^ 7^ 0. 

Theorem 3.2. Suppose that pi, . . . ,pn is a gci and pn = xf^ — c„a;™_]^, with £,m 
coprime positive integers. Letpi, . . . ,Pn-i be the binomial system obtained through 
parametric reduction. Then D = D and, for any L <Z [n], 

(3.2) ^L=/i£. 
Consequently, d = d as well. 

Proof. Let c = (ci,...,c„) G (fc*)" be such that Jc is a complete intersection. 
We may assume without loss of generality that c„ = 1. Let c = (ci, . . . , c„_i) 
and denote by Jc the ideal generated by pi(c;u), . . . ,p„_i(c;u) in the ring k[u\. 
Given any A = (Ai, . . . , A„_i) G ^{Ji) C fc"~^, let us denote by A the point 
(Ai, . . . , A„_2, A^_i, A™_j) e 'V{J'c) C This assignment A 1-^ A defines a bijec- 
tion between N{Jc) and V(j7c) since m are coprime, and so D = D. To show that 
d = d it suffices to prove that at the level of local rings 

(3.3) dim-k{R(g)kk)x/{Jc)x = dim^(i? ®fe fc)^/(j£)3^. 

We will denote by Ai the localization of k[ui, . . . ,u„-i] at A and by A2 the 
localization of k[ui, . . . , Mn-2, "l-i, at A. Let (Js)^^ be the ideal generated 

by pi{c;u), . . . ,pn~i{c;u) in A2 so that {Jij^ — ^1 ' (^7^)3^• Again, since m and £ 
are coprime it is clear that 

diiaj^{R®kk)\/{.Jc)\ = Aim-^A2/{Ji)^^. 
Thus, the result will follow if we show that 

(3.4) dxTn-^A^/iJ,)-^ ^ dnT,-^A2/{Jc)i 

The following proof of (|3.4|l was suggested to us by Mircea Mustafa. 

We recall from j21[ §14] the following notion of multiplicity: Let (i?, m) be a d- 
dimensional Noetherian local ring, M a finite i?-module and q an m-primary ideal. 
The multiplicity of M with respect to q equals 

(3.5) e(q,M) = lim length(Af/q™+iM) 

m — yoo rn 

Since both Ai and A2 are Cohen-Macaulay rings of dimension n ~ 1 and 

Pi{c;u), . . . ,_p„_i(c; u) 

define a regular sequence in A2, hence in Ai as well, it follows from '21', Theo- 
rem 14.11] that 

dlm-^A2/{Jc)~^ - t{{Jc)i.A2) and dim^ ^1/(^5)3^ = e(( Jg)^,, Ai). 
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On the other hand, Ai may be considered as a A2-niodule and it is clear from H3.5|l 
that 

Finally, [HI Theorem 14.8] gives that 

e((Jc)A,^i) =rankA,Ai •e((j£)3,,A2) = e(( Jg)^^, ^2), 

since the assumption that m and i are coprime implies that the two domains Ai , A2 
have the same fraction field and so rank^^Ai = 1. This proves H3.4|l . □ 

3.2. Irreducible Systems. 

Definition 3.3. A binomial system pi, . . . ,p„ is said to be irreducible if it is in 
normal form and it is not possible to reorder it so as to find a proper index subset 
/ C [n] such that for every i G I the binomial pi depends only on the variables 

xj,j e /. 

Recalling that a system in normal form is a gci and that G V(J'), we easily 
have: 

Lemma 3.4. Let pi, . . . ,Pn he an irreducible system as in and let c € (fc*)" be 
such that Jc is a complete intersection. Then if a Wc, either a = or a G (fc*)". 

Proof. Given a e Vc, let / = {i G [n] : 7^ 0}. If i G / then, since pi, . . . ,p„ is in 
normal form, 

Pi(c; x) = xl' - CiX^' ; ri > 0, (3i ^ 0, 
and, since a; ^ 0, we must have supp(/3i) C / for all i £ I. This contradicts the 
irreducibility of pi, . . . ,p„ unless / = [n] or / = 0. □ 

The following theorem identifies d and /i for irreducible systems. Recall that 
6 = I det B\ is the cardinality of Vc n (A:*)". Our arguments are built on the proof 
of a result of Vinberg (cf. (18i Theorem 4.3]). 

Theorem 3.5. Given an irreducible system 

Pi{c;x) = x\' - CiX^^' , i = l,...,n, 
where ri > 0, f3i e W\ f3i ^ 0, then: 

• If all principal minors of B are positive 

d ~ ri ■ ■ ■ rn ; fi = d — S. 

Such a system will be called a global irreducible system. 

• Otherwise, /i = ri ■ • ■ r„ and d = fi + S. In this case we say that the system 
is local. 

Proof. Let us fix throughout coefficients c G (fc*)" such that J7c is a complete 
intersection. Since the system is in normal form, the entries of B are ha — ri — {[3i)i 
and bij = —{Pi)j,i ^ j. Hence, its off-diagonal terms are non-positive. Moreover, 
the irreducibility of the system implies that B is indecomposable in the sense of 
[18|. In fact, the irreducibility of the system implies a stronger condition, namely 
[181 Lemma 4.3]: Suppose u G M" is a vector with non-negative entries and that 
B-u> in the sense that all its entries are non- negative as well. Then either u = 0, 
or u > 0, i.e., all its entries are strictly positive. Indeed, let I = {i £ [n] : Ui — 0, 
then for any i G I, {B ■ u)i < and equality occurs if and only if bij — for all 
j ^ I. Hence, by irreducibility we must have / = [n] or / = 0. 
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Given that |18l Lemma 4.3] holds in our case, we can apply Theorem 4.3 in ^H] 
and conclude that three cases are possible: 

• There exists w G Q" all of whose entries are positive such that B ■ w > 0. 

• There exists w G Q", all of whose entries are positive such that B ■ w < 0. 

• rank(_B) — n — 1 and there exists w G Q" all of whose entries are positive 
such that B ■ w — 0. 

According to T, Theorem 2.3], the first condition is equivalent to the statement 
that all principal minors of B are positive which implies, in particular, that all the 
diagonal entries of B are strictly positive. These are the so-called M-matrices of 
Moreover, if we consider a term order in k[xi, . . . ,Xn] that refines the weight 
order defined by w, the term x^' will be the leading term in pi{c;x), and hence 
Pi(c; x), . . . ,p„(c, x) is a Grobner basis. It then follows '3', §5.3, Proposition 4] that 
d = ri ■ ■ ■ Vn and, by Lemma 13.41 fi = d — S. 

In the second case we can similarly define a local order (cf. JA^) for which the 
leading term of (c; x) is x^^ . Hence pi{c; x) , . . . , pn{c, x) is a standard basis in the 
local quotient ring at the origin and, consequently, /i = ri ■ • ■ r„ and d = /i + 5. 
We note that this is valid even if det B = since, in that case, J7c is a complete 
intersection if and only if Vc = {0}. 

In the third case, the binomials pi{c]x) are weighted homogeneous relative to 
the weight w and therefore /i = ri • • • r„ and d = /i + <5 since, again, Vc consists of 
only the origin. Thus this case behaves as the previous one and we will also refer 
to it as a local case. □ 

Remark 3.6. We note that if n = 1, the system p = x°' — cx^ ,a ^ (3, will be local 
if a < /? and global if a > /3. 

3.3. The General Case. We consider now general gci systems in normal form. 
Throughout this subsection we will, again, fix coeSicients c G (fc*)" so that Jc is 
a complete intersection. For economy of notation we will denote simply by pi the 
corresponding binomials in fc[xi, . . . , a;„]. If the system pi, . . . ,p„ is not irreducible, 
then, as Lemma l3 . 81 shows . it is possible to choose an increasing sequence 

(3.6) = 1^0 < < ■ ■ ■ < = n 

so that if la — {I'a-i + 1, • • ■ , i^a}, then the following holds: 

• For i e /a, Pi e k[xj;j G /i U • • • U /a]. 

• The system pi :— pi(l, . . . , 1, . . . , x,^^), i G la, is irreducible. 

Definition 3.7. A system of this form will be said to be in triangular form relative 
to the blocks /i, . . . , /s. Given a reducible system in triangular form, we will refer 
to the system {pi, i £ la} as the restriction of pi, . . . ,p„ to la and denote it, for 
short, by p". 

Lemma 3.8. Any system of n binomials pi, . . . ,p„ in normal form f2.9\) can be 

put in triangular form in time 0(n^). 

Proof. Consider the ocurrence matrix N: this is a 0-1 matrix with 7^ if 
and only if i j and pi depends on Xj (i.e., if pi = x^' — CiX^^ with Pij ^ 
0). This is a standard construction, first used by Steward |23, for the analysis 
of the structure of large systems of equations. Note that, because the system 
is in normal form, putting pi , . . . , p„ in triangular form corresponds precisely to 
finding a permutation matrix P such that *PNP is block lower triangular, with 
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the irreducible subsystems of pi , . . . , p„ corresponding to the irreducible diagonal 
square blocks along the diagonal of *^PNP. 

Tarjan's algorithm '31? to search for the strongly connected components of the 
directed graph associated to N provides an efhcient method for finding such per- 
mutation matrix P [2| I23| ; it runs in time linear in the number of vertices plus the 
number of edges of the graph. □ 

Given a system in normal form and triangular relative to Ji , . . . , Jg , let 5a = 
I dot Ba I , where Ba is the matrix associated with the system p° and 



Pa 



n 

3 1^1 a 



We also denote by /Iq the multiplicity of p° at and by da the total number of 
solutions of counted with multiplicity. 

For a triangular system pi , . . . , p„ , its associated matrix is block lower-triangular: 



(3.7) 



B 



( 

C21 





B2 



\ 





\ Csi Cs2 ■ ■ ■ Bg J 

The number of solutions of the system pi , . . . , p„ and the patterns of possible 
zero coordinates of the solutions are best described in terms of the directed acyclic 
graph G with s vertices labeled {1, . . . , s} and an arrow from node a to node b if 
and only if the rectangular submatrix Cba is not identically zero. We recall that a 
vertex is called a source if it is not the head of any arrow. The subset of sources 
of the vertex set of a subgraph H oi G will be denoted by S{H). 

Remark 3.9. We can think of G as a weighted graph, where each vertex a G [s] 
comes with the weights Sa^Pa^Pa (or Sa^da, ^la)- Equivalently, we can think that 
the information at each node is coded by the weights 5a^ Pa plus an additional label 
global or local according to where Ba is global or local, which prescribes the relation 
among 5a, Pa and p.a (or 5a, da and pa)- 

Theorem 3.10. The multiplicity p of Jc at the origin equals 



n 

yaeG\S(G) 




(3.8) 



Proof. We will prove formula 13. 81 bv induction in the number s of blocks. If s = 1, 
the system is irreducible and {1} G S{G) so the formula holds. Consider s > 1 and 
assume that the result is true for systems with s — 1 blocks. Let B be as in H3.7|l . 
set n' :— Vg-i, where Vg-i is as in 1)3. 6|l . and consider the ideal := (pi, . . . ,Pn') , 
in the polynomial ring in the first n' variables. Clearly, pi, . . . ,p„' is in normal and 
triangular form. Let G' be the corresponding graph; it is obtained by erasing from 
G the vertex s and all edges ending at s. By inductive hypothesis, we have that 
the multiplicity p' of J^. at 0' equals 



(3.9) 



n 

ya(^G'\S{G') 
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\ 


B' 





C 


Bs J 



The matrix B has the form 

/ 

(3.10) B = 

V 

If the rectangular matrix C is identically zero, then the last n — n' polynomials 
depend only on the last n — n' variables, and we have that 

= • Ms, 

as wanted, since in this case S{G) = S{G') U {s}. 

On the other hand, if C is not zero, it is possible to find a positive weight vector 
w such that the initial monomial in-w{Pj) — Xj\ for all n' < j < n. Indeed, set 
Jo = [n'] and, for I > 1 define 

Je := i^ke [n]\ (^[j j}j : Ji-i D supp(/3fc) ^ 0| . 

Note that C ^ implies that Ji is non empty. Also, the assumption that Bg is 
irreducible guarantees that there exists L < n — n' such that — Ui<£<l J^- 

Now, choose = 1 for /c e Jl. Then assuming that the weights for the variables 
k S Ja, £ < a < L — 1, have been chosen so that in^^{pj) — a;^^ for all j S J;,, 
b > we may choose positive weights Wk for k G Je-i that are sufficiently large 

so that in-w{pj) — x^f for all j S Ji as well. 

Consider now any local order -< in k[xi, . . . ,a;„] refining the weight —w. Let 
{qi, ...,qt} be a standard basis for the ideal J'^ with respect to the local order 
induced by -< in k[xi, . . . , Xn']- Then, {gi, . . . , qt,Pn'+i, ■ ■ ■ ,Pn} is a standard ba- 
sis for J^c relative to -< since, for every i — 1, . . . ,t, the leading monomials of the 
polynomial qi is coprime with those of the pj, n' < j < n, and, therefore, the weak 
normal form of the corresponding S'-polynomial is |14| . The corresponding initial 
ideal L^{Jc) will be generated by some monomials in the first n' variables (gener- 
ating the initial ideal L^i{J'^j) and the pure powers x'j^ for all j > n' . Therefore, 
the multiplicity /i of J^c at equals: 

dinifc (k[xi . . . ,a;„]/J"c)o = dinife (k[xi . . . ,Xn]/L^{Jc))^ = 

dinifc {k[xi . . .,Xn']/L^,{X))^ ■ dimj, {k[xn'+i ■ ■ ■ , a;„]/(a;^7+Y . ..<"))(, = 

dinifc {k[xi . . .,Xn']/X)^ ■ Ps- 

In this case s ^ S{G), and so S{G) = S{G'). Since the dimension of the local 
quotient by J'^ at the origin equals H3.9|l . we get that 

M = m' ■ Ps n ^'^ n ^'^ 

\aG\S(G) / \b£S{G) 

as wanted. □ 
Remark 3.11. Using Theorem 13.51 we can translate H3.8|) as 

(3.11) f= ili'^-] in . 
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where Gi is the set of nodes of G corresponding to the global, non-sources of G 
and G2 is its complement. 

We will also need the following terminology. 

Definition 3.12. A vertex h of (the directed acyclic graph) G is said to be a 
descendant (respectively, a direct descendant) of the vertex a if there is a directed 
path (respectively, a directed edge) from a to &. A (directed) subgraph iJ of G is 
said to be full if, for any of its vertices j, all its descendants and all the directed 
paths starting from j also belong to H . The collection of full subgraphs of G will 
be denoted by T{G). 

The empty subgraph is full and even if G is connected, a full subgraph H may 
be disconnected. Note also that a full subgraph is completely determined by its 
sources. 

The following result refines the description given in Remark [2.7l of subsets L C [n] 
with /ii 7^ . 

Proposition 3.13. Let pi, . . . ,pn he a binomial complete intersection in normal 
and triangular form and L C [n]. Then fi^ ^ unless there exists a full subgraph 
H of G such that 

(3.12) n ^ 

and L coincides with the union of all the indices belonging to blocks that are vertices 
ofH. 

Proof. With the above notations, let A = (Ai, . . . , A„) G V(j7'c) and L = L{X) = 
{i G [n] : Xi — 0}. Set = {a S G : /a n L ^ 0}. U a £ H then we may argue as 
in Lemma [3.41 to conclude that la C L. Suppose now that a £ H and that (a, b) is 
an edge in G. Since Cba 7^ 0, there exists i & la and j G h such that i G supp(/3j) 
and, consequently, Xj = 0, i.e., j G IbC] L, and b G H. This shows that _ff is a full 
subgraph of H. The need for condition (|3.12|) was already noted in Remark 1^71 □ 

With notation as in Prop. lTTSl given a full subgraph H C G, we will denote by 
L{H) the set of indices belonging to blocks associated with vertices of H. 

Proposition 3.14. Given a full subgraph H ofG, the number D]^(^h) of points in 
V(j7c) n k2(H) counted without multiplicity equals 



(3.13) 




while the number fJ.L{H) of points in V(j7c) ^^^l(h) counted with multiplicity equals 



(3.14) - 




Proof. The first assertion follows easily from Proposition 13.131 In order to prove 
(EH, let A e V(Jc) n write A = (A^^), . . . , A^")) with A^'') G fcl^-l for all 

a G [s]. Since _ff is a full subgraph, there are no edges starting at a node in H and 
ending at a node outside of H; i.e., Cba ~ for all a G and b ^ H. Therefore, it 
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is possible to relabel the variables and the binomials pi , . . . , p„ so that the system 
remains in normal form and satisfies that a < b for all a ^ iJ and b Cz H. Thus, 
we may assume without loss of generality that H = {t + 1, . . . , s} and therefore 
A= (A(i),...,A(*\0,...,0) with A^'^) e (fc*)!^"! for a = Equivalently, 

A = (A',0) e (fc*)"' X (fc)"-"' ; n' i^t • 

We let x' stand for the first n' variables xi , . . . , Xn' and x" for the remaining n^n' 
variables. Then 

J'. := {pi, . . . ,p„'} C k[x'] 

and A' is a simple zero of J^. Hence pi, ■ ■ ■ ,Pn' define the maximal ideal in the 
local ring {k[x'])^,. We then have: 

fj,\ := dim^ (fc[a;]/J"c);^ 

= dimfe (k[x]/{xi ~ Xi, . . . ,Xn' ~ Xn',Pn' + l, ■ • ■ ,Pn))^ 

- dimfe {k[x"] I (p„,+i ( A' , x" ) , . . . , p„ ( A' , x" ) ) ) 

- dimfe (fc[a;"]/(Pn'+i(l, ■ • • , . • . ,P„(1, . ■ . , 1, x")))^ . 

So, pu\ equals the multiplicity at the origin G fc"^" of the system {p„'+i, . . . 
Formula (|3.14(l now follows from Theorem 13.101 and the fact that the system 
Pi, . . . ,p„' has 5i • ■ • (5t simple solutions in (fc*)" . □ 

The following explicit formulas for d and D follow by adding H3.13|l and (|3.14|l 
over all full subgraphs of G. 

Theorem 3.15. Suppose that pi, . . . ,Pn are in normal, triangular form. For 
generic parameters c £ (fc*)", the total number of solutions of the system pi(c; x) = 
■ ■ ■ — Pn{c; x) = 0, counted without multiplicity, equals 

(3.15) D= J2 (llM, 

and the total number of solutions counted with multiplicity equals 



(3.16) ^= E (iim ( n p" 

Her{G) \a<^H ) \beH\S{H) 

We end this section with a recursive formula to compute d. In order to state the 
following proposition we define, for 1 < r < s, the binomial system g^*"^: 

Pi{l, . . . , 1 

Note that the matrix associated with q^''-' is: 




(3.17) Sf'') 



/Br ... \ 

C(r+l)r Br+1 ... 



\ Csr Cs{r+1) ■ ■ ■ Bs J 

Clearly if pi, . . . ,p„ is in normal, triangular form, so is q^^\ We denote by the 
number of solutions in counted with multiplicity, of the system q'^^\ 



18 EDUARDO CATTANI AND ALICIA DICKENSTEIN 

Proposition 3.16. Fr is a polynomial function of {5a, ^ia, Pa] a = r, ...,s}. It 

may be computed recursively as: 

Fs ^ ds ^ 5s + Us 



(3.18) Fr — 5r ■ Fr+1 + ' i^r+1 1 i56=0,A»6=p6 j 

where b runs over all indices in {r + 1, . . . , s} such that Cbr 7^ 0. 

Proof. We may assume without loss of generality that r = 1 < s. Let G be the 
graph of B and G*^^^ the subgraph of G associated to the submatrix B^^^ defined 
by <TT7ll . 

Any full subgraph H g J^(G^^-') may be thought of as a full subgraph in G. We 
denote by T' C J-{G) the collection of such subgraphs. Clearly T' consists of all 
full subgraphs of G not containing the vertex 1. Let T" denote the complement 
of JT' in J-{G). Removing the vertex 1 from a subgraph H e IF" defines a full 
subgraph i/*^^) of G^^^ with the property that no direct descendant of 1 in G may 
be in G^^^ \ H^'^\ Let us denote by F"{G^'^^) the collection of such full subgraphs 
of G'^' . We can write 

(3.19) ^1 ^ X! l^L{H}+ l^L{H}- 

Since, for H G !F', 1 ^ iJ, in view of (|3.14|) . the first sum may be computed as: 

(3.20) ^ ^LiH) ^ 5i ^ fiL{H) = 5iF2, 

Hey^' HG.F(G(2)) 

since S{H) is the same whether we view H as a. subgraph of G or of G^^-*. 

Thus, in order to complete the proof we need to show that the second sum in 
(|3.19|l equals 

Ml • F2\st=o,iib=Pb^ 

where b runs over all vertices in G'^^ that are direct descendants of 1 in G. We 
note first of all, that setting 5^ — for all direct descendants 6 of 1 has the effect 
of restricting the sum in (pTT^ to J^"{G^'^^). Moreover, given H e F" , let i/^^) 
denote the full subgraph of G^^^ obtained by removing the vertex 1 from H. Then 
S{H'-^'>) consists of S{H) n G^^^ together with all the direct descendants of 1 in H. 
This change may be accomplished by replacing fXb by pt whenever b & H^"^^ is a 
direct descendant of 1 in H. Since 1 G S{H) for all H € F", we obtain the desired 
equality. □ 

Example 3.17. We return to Example 12.131 We recall that the reduced system 
Pi,...,P6 is: 

pi^ul- u\u2] P2^ul- ulu2; 
P3^ul~ulul; p4 ^ ul - ulu2ulue; 
P5 = ul-ulu5; pe = ul ~ uluiU^UQ. 
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and, therefore, its associated matrix is 

\ 





1 



2 / 

Therefore, the system is in normal, triangular form with blocks relative to the index 
sets /i = {1, 2}, /2 = {3}, and I3 = {4, 5, 6}. The block Bi is global, while B2 and 
B3 are local. The graph G has 3 vertices {1, 2, 3} and arrows from 1 to 2 and 1 to 
3. Hence S{G) = {!}. The weights are: 

<5i = 4, ^2 = 1, <53 = 5, pi = 35, p2 = 2, p3 = 18 , 

and, taking into account the local/global label, we get ^1 = 31, /12 = 2, /is = 18. 

We may now apply 1)3.8(1 to compute the multiplicity /i of (pi, . . . ,pg) at the 
origin: 

A = Ml • P2 • P3 = 1116 . 
In order to compute d we use the inductive procedure of Proposition l3.16l Since 
the subgraph with vertices {2,3} is disconnected we have: 

F2 = {S2 + M2) • {h + Ma)- 

Hence, Fi — Si ■ {62 + M2) • (1^3 + Ma) + Mi ■ P2 • Pa- This gives d = 1392. We note 
that this is far from the Bezout bound of 43740. 

Using Lemma ITTI and Theorem 13. 21 we see that the total number of solutions for 
the original system pi, . . . ,p8 are given hy d = 2d and m = 2/1. This values may be 
easily verified using a computer algebra system such as Singular . 

Finally, we note that G has five full subgraphs with vertex sets: {1,2, 3}, {2, 3}, 
{2}, {3}, and 0. This means that there are five index sets L C [6], such that /i^ 7^ 0. 
They are Li = [6], L2 = {3,4,5,6}, L3 = { 3}, L 4 = {4,5,6} and L5 = 0- The 
corresponding multiplicities are according to 1(3.14(1 : 

Mz, = M = 1116, Ml. = 144, ML3 = 40, Mz, = 72, Ml, - ^ = 20. 
Moreover, the total number of solutions counted without multiplicity is given by: 
D ^ 61 + 61 ■ d:i + Si ■ 62 + Si ■ 52 ■ S3 = 48 

This information may be lifted to the original system using the bijection L ^ L 
discussed before Theorem 13. 21 We get that ml = except for the following subsets 

Li = [8], L2 = {4, 5, 6, 7, 8}, L3 = {4}, L4 = {5, 6, 7, 8}, L5 = 0- 
Once again, fi^. = 2 ^ij^ . 

4. Counting complexity 

In this section we will study the counting complexity, in the sense of (31( . of com- 
puting the numerical invariants d, D, S, m, ^ind ml associated with a gci pi, . . . ,Pn- 

We have already proved that we may decide in polynomial time if pi , . . . , p„ is 
a gci and that the property of being a complete intersection is independent of the 
coefficients if det B 0. Moreover, if pi, . . . is a gci we may also transform it 
into normal and triangular form in quadratic time. Also, since a system with generic 
exponents is irreducible and satisfies detS ^ 0, we may compute its invariants in 
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time polynomial in n for any choice of coefficients by Theorem 13. 51 In the general 
case, we may compute S, /i, and /i^, for a particular choice of L, directly from the 
invariants 5q, pa, and /i^ associated with the diagonal blocks of the system. Thus, 
(5, /i, and fj,L may be computed in polynomial time as well. 

However, we will show below in Theorem 14.31 that the computation of d or Z? 
is a #P-complete problem, and therefore it is at least as hard as an NP-complete 
problem |31j . In order to do this we begin by reversing the relationship between 
binomial systems and weighted acyclic directed graphs. We recall that to a bi- 
nomial system , . . . , p„ in normal and triangular form we associate an acyclic 
directed graph G whose vertices {1, . . . , s} correspond to the diagonal blocks of the 
associated matrix B and that each vertex has weights 6a, Pa, o, [s], plus a label 
"local" or "global" . In the first case we set Pa — Pa, while in the global case we set 
Pa = Pa — 6a- In any case da — Sa + Pa- The proof of the following proposition is 
straightforward. 

Proposition 4.1. Let G = {V,E), V = [s], he an acyclic directed graph, with 
weights 5a, Pa G ^>o o.nd labels local/global attached to each vertex. Let pa and da 
he defined as above. Then, the system of binomials defined by 



for all local vertices a, has as weighted graph {G,Sa, Pa, Pa)- 

Remark 4.2. The total number of solutions d and D of the system in Proposi- 
tion O are given by (|3.16|l and (|3.15() , for generic parameters Cq . For any order on 
the set of vertices of G such that i < j if there is a path from node i to node j (i.e., 
for any linear extension of G), it is clear that the corresponding matrix B of the 
system will be lower triangular, with diagonal entries ±(c?a — pa). Thus, whenever 
da 7^ Pa, we have that det(_B) ^ and we may simply choose = 1 for all a G [s]. 

Note also that if a is a source of G, then we get Pa — xf" — CaX^^'' in the global 
case, and Pa = x^" — CaX^" in the local case. This is compatible with Remark 13.61 

In the particular case when all vertices {1, . . . , s} of a directed acyclic graph G 
are local, and their weights are (5a = 1, Pa = 1, for all a G [s], the binomial system 
defined in ProDOsition l4. II takes a very simple form: 



We will refer to this system as the standard binomial system associated with G. 

Theorem 4.3. Computing d and D for binomial complete intersections pi, . . . ,p„ 
in normal, triangular form are ^P-complete problems- 
Proof By Theorem 13. 151 the problems of computing d and D are in the complexity 
class #P. We will show that computing these invariants gives, for special binomial 




for all global vertices a, and 




(4.1) 
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systems, the number of independent subsets of a bipartite graph G. Since, by |24j . 
this is known to be a #P-complete problem the result will follow. 

Let G be a bipartite graph with vertices { 1 , . . . , s} . Let pi , . . . , be the standard 
binomial system of G as in H4.1|) . Then, for each full subgraph H C G we have, by 
(|3.14|l . that Hl{h) — 1- Hence, according to H3.16|l and H3.15() . both d and D are 
equal to the number of full subgraphs of G. But, as has been noted earlier, a full 
subgraph is completely determined by its sources and, for a bipartite graph G, a 
subset of vertices is the set of sources of a full subgraph H if and only if it is an 
independent subset of G. Thus, d and D agree with the number of independent 
subsets of G. □ 

Recall that a directed acyclic graph G = {V, E) is called transitive if there is 
an edge (a, 6) G E each time that there is a directed path from a to 6, Transitive 
directed acyclic graphs are in correspondence with partial orders -< on V , where 
a < b \i and only if (a, 6) G E. Given a partial order on a subset A of F is 
called an antichain if given ai,a2 £ A, neither ai ^02, nor 02 -< ai. It is shown 
in |24| that counting the number of antichains in posets is a #P-complete problem 
and, hence, #P-hard. Given any directed acyclic graph G — (y,i?), it is possible 
to compute its transitive closure G+ = (V, E^), in time 0(|V^p) by the well known 
Floyd- Warshall's algorithm. It follows from (|3.16() and (|3.15|) that d and D are the 
same for the standard binomial systems associated with G and with G+. 

Proposition 4.4. The number of (simple) solutions of the standard system j^. J| ) 
associated with a directed acyclic graph G equals the number of antichains in the 
associated partial order. 

Proof. As in the proof of Theorem 14.31 for the standard binomial system of G 
we have d = D and this number agrees with the number of full subgraphs of G. 
These subgraphs are determined by their sources, which correspond exactly to the 
antichains in the associated partial order on V . □ 

Although, as the previous results show, the problem of computing the total 
number of solutions for a general binomial system in normal and triangular form is 
7^_P-hard, there are classes of binomial systems whose invariants may be computed 
in polynomial time. For example, if the graph is totally disconnected then d = 
di - ■ ■ ds = Y\l^i{5i + fii). At the other extreme if G is a (complete) directed graph 
with vertices {1, . . . , s} and {b, a) is an edge of G for all a.b E [s] with a < b, then 
it is easy to see that there are only s + 1 full subgraphs of G and, consequently, the 
sums in H3.15|l and (j3.16(l . consist of s + 1 terms. 

Even if the number of full subgraphs is exponential in s and G has few connected 
components, a bound on the number of local blocks guarantees that d can be 
computed in polynomial time in ri. For instance, if all blocks are global, then B is 
an Af- matrix and pi, . . . ,p„ is a Grobner basis for a positive weight order, and so 
d — pi ■ ■ ■ ps. We end with the following "positive" complexity result. 

Proposition 4.5. Let N G Z>o. Assume pi, . . . ,p„ is in normal and triangular 
form with s blocks of which at most N are local. Then, there is a formula to compute 
the total multiplicity d with at most 2^ summands, each involving s products. Thus, 
if the number of local blocks of a binomial system in normal and triangular form 
is bounded independently of n, the number of affine solutions of the system can be 
computed in time polynomial in n. 
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Proof. Recall the notation in Proposition 13. We may write the polynomial for- 
mula Fr{{5a, Pa, ^J■a)^ a- G [s] ) for the computation of the total number of solutions 
of the system purely in terms of 6a and pa by keeping track of the local/global 
character of each vertex and replacing pa by Pa if a, is local and by Pa — Sa in the 
case of a global vertex. We call Fr{{Sa, Pa),^ G [s]) the polynomial obtained after 
these substitutions. Then, for a global vertex r, the recursion H3.18|l becomes 

(4.2) Fr^6r- Fr+1 + [pr - 5r) ■ /'^+i|5„=0, 

where a runs over all direct descendants of r. Let us write f'r+i = -F^+i + 
where i^^+i consists of all summands containing a factor 5a with a a direct descen- 
dant of 1. Hence, F^-i-i vanishes when we set such 5a — Q and H4.2|l becomes: 

Fr^5r- + F'r+l) + {pr - ^r) ' = ' + Pr ' -F'r'+l 

and, consequently, the total number of summands does not change when adding a 
global vertex. 

On the other hand, if is local then (|3.18|l becomes 

Fr — 5r ■ Fr+l + Pr ' -Pr+l|(5„=0 

and the number of summands is, at worst, doubled. 

It follows that when A'' is bounded independently of the number n of variables, 
d can be computed by adding a constant number of summands. Each of these sum- 
mands has s < n products of factors involving the computation of determinants of 
the square diagonal blocks of the associated matrix B or products of the exponents 



5. Applications 

In this section we will briefly discuss some of the problems that led us to the 
study of systems of n binomials in n variables. 

An important subfamily of binomial ideals is given by the toric ideals associated 
to configurations A — {ai, . . . , a„i} C Z*"' of integral points spanning Z*^: 

Ia = (a;'" -x'' ; A-{u-v) = 0), 

where u,v G N™. In particular, beginning with the work of Herzog and De- 
lorme [21 the question of classifying complete intersection toric ideals (and the 
corresponding semigroup algebras) has been extensively studied by many authors 
m 01 ^2 1121 IIHI |2E| ■ A key step in many of these works is the study of the ideal 
generated by binomials a;"' — a;"' associated with a Z-basis of the kernel of A. 
More generally, given Q-linearly independent elements , . . . , j/,. G , consider 
the associated lattice basis ideal J C A:[a;i, . . . , Xm], generated by the binomials 

bj ^ x''^ ~ x"^ ; j = l,...,r, 

where Vj = uj — Vj, and Uj,Vj G N™ have disjoint support. Let C C Z'" denote 
the lattice spanned by i^i, . . . ,i>r and let Ic ■= (x'^ — x"" : u — v € C) he the 
corresponding lattice ideal. We assume that these ideals are homogeneous, i.e., 
wi + ■ ■ ■ + Wm = 0, for every w & C. 

The ideal Ic is prime if and only if the lattice C is saturated. If C is not 
saturated, then Ic has g radical primary components, where g is the index of C in 
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its saturation. Moreover, all these components have the same degree, equal to the 
degree dc of the associated toric variety ^O] • 

We can apply Theorem 13.151 to compute the multiplicity and geometric degree 
[21 of the primary components of J. This may be used to describe the holonomic 
rank of Horn systems of hypergeometric partial differential equations and to study 
sparse discriminants, generalizing the codimension-two case |H1[Z|. 

A straightforward extension of the results of 17 to non-saturated lattices gives 
the following description of all primary components q of J. Let K C {1, . • . 
and Z(K) C {l,...,r} as in (|2.3|l . Assume that n := \Z{K)\ = \K\ and for aU 
] i ZiK) 

supp(uj) n K = supp(wj) n x = 0. 

Let p' be a primary component of the lattice ideal Ic' associated to the sublattice 
of Z"-" spanned by Vj, j ^ Z{K). Then, the ideal 

q = p' + {h,i ^Z{K)) 

is a primary component of J with associated prime 

p = p' + {xk,keK). 

Note that for i^T = we recover the components of Ic- 

In order to describe the multiplicity and geometric degree of a component q, 
let us assume that K — Z{K) — {1, . . . ,n} and for any w G Z™, denote tt{w) = 
{wi, . . . , Wn)- Let aj — 7r(uj), f3j — ■n(vj) and set 

Pj{c;x) = x"^ — CjX^^ ,Cj e k* . 

Since J is a complete intersection, pi , . . . , p„ is a gci. Let ^ denote the multiplicity 
at the origin. Fix coefficients c G (fc*)" such that Jc is a complete intersection. 
Since 

^ = length {k[xi, . . . ,x„]/Jc)o = length (fc[a;i, . . . ,a;™]/J)p , 
and the degree of p equals that of p', we have 

Proposition 5.1. With notation as above, the multiplicity of q equals fi and the 
geometric degree of q equals dc> ■ fJ-- 

As a second application, consider a system of constant coefficient partial differ- 
ential equations defined by n operators of the form 

(5.1) a,9"^-6,a^- j = l,...,n, 

where aj,bj S k* , aj,f3j £ N", aj ^ /3j. Assume moreover that the ideal J in 
k[xi, . . . ,Xn] generated by the binomials ajx"'^ — bjX^^ is zero-dimensional. As 
before, let fi^ the number of points in V( J) H fc£ counted with multiplicity. From 
[291 Chapter 10], we have the following characterization. 

Proposition 5.2. Let L C {l,...,7i}. The dimension of the space of solutions 
to which depend polynomially on the variables xi, £ £ L, and exponentially 

on the remaining variables Xj, j ^ L, equals fi^. 

These dimensions can then be computed using the results in Section O particu- 
larly formula (|3.14() . 
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